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Over the years, a variety of the group-theoretic methods for the study of 
differential equations have been devised [1,2]. It is scarcely to belabor their 
importance in the investigation of integrable systems. Let us only recall the 
methods for constructing group-invariant solutions. Nevertheless, there is still 
a general lack of approaches discussing behaviour of nonintegrable equations. 
Recently, an interesting approach was proposed by Oknihski [3] relying on 
the study of nonlinear recurrences by asscociating them with a discrete-time 
evolution on a Lie group. More precisely, within the Oknihski approach the use 
is made of the Shimizu-Leutbecher sequence [4,5] which satisfies the following 
group recurrence: 

9k+i = gkhgk 1 , (1.1) 

where g^s and h are elements of a Lie group G. On taking a concrete re- 
alization of a Lie group G we arrive at the system of nonlinear difference 
equations which obey the group parameters. For example in the case with the 
577(2) group one obtains the logistic equation satisfied by one of the coordi- 
nates of the system implied by (1.1). By choosing G = E{2) and G = SU{2) 
for the Schimizu-Leutbecher sequence (1.1) the method was applied to the 
study of random walk on the plane [3] and discrete symmetries in chaos-order 
transitions [6], respectively. The question naturally arises as to whether such 
approach can be generalized to involve the case with continuous-time equa- 
tions. An attempt in this direction was made by Oknihski [7] who studied 
the equation with delayed argument arising from (1.1) by setting k to be a 
continuous parameter. 

In this paper we introduce an abstract Newton-like equation on a general Lie 
algebra which can be regarded as a continuous-time algebraic version of (1.1), 
such that the orbits of the Lie-group action are attractors. The correspond- 
ing group parameters fixing a concrete algebra obey the system of nonlinear 
second-order differential equations having a limit set coinciding with the orbit. 
In section 2 we introduce the abstract Newton-like equation on a Lie algebra 
and study its asymptotic behaviour. We show that orbits of the Lie-group 
action are attracting set. Section 3 discusses periodic solutions to the abstract 
Newton-like equation. Namely, we establish some criteria of the existence of 
such solutions. In section 4 we illustrate the theory by an example of the 
SU(2) group. We first discuss the nonlinear dynamical system generated by 
the corresponding Newton-like equation on the su(2) algebra and its asymp- 
totic counterpart. We then study the dynamics on the orbit i.e. the sphere 
S 2 which is the limit set for the trajectories of the discussed nonlinear sys- 
tems. Such the dynamics is demonstrated to correspond to a generalization of 
the spherical pendulum involving the presence of the nonpotential forces. We 
identify the periodic solutions and show that the considered dynamical system 
on the sphere has the global bifurcation of a limit cycle from an equilibrium 
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mentioned above we conclude that this system has the global Hopf bifurcation. 



2 The Newton-like equation 



In this section we introduce the abstract Newton-like equation on a general Lie 
algebra. By choosing the concrete algebra this equation generates the nonlinear 
dynamical system satisfied by the time-dependent group parameters. Consider 
the following second-order differential equation on a Lie algebra $j: 

pX + uX + pX + aY = e iX Ye- iX , X(0)=X o , X(0) = X , (2.1) 

where X(t): R — > is a curve in g, Y e is a fixed element, p > 0, v, p, a e R 
and dot designates differentiation with respect to time. 

Demanding that (2.1) admits the solution on the orbit 

TrX 2 = const ^ 0, (2.2) 

assuming that p, v — const and p, a are independent of Y, and using the 
differential consequence of (2.2) such that 

TrXX = -TrX 2 , (2.3) 

we find 

TrX 2 

a=1, P = ^TX^- (2 ' 4) 

Inserting (2.4) into (2.1), rescaling t — > y/pZt and setting = (5 we finally 

obtain the following Newton-like equation: 
TrX 2 

X + pX + —X = e lX Ye-* x -Y, X(0) = X , X(0) = X . (2.5) 

We now discuss the asymptotic behaviour of (2.5). An immediate consequence 
of (2.5) is 

TrXX + (3TrXX + TrX 2 = 0. (2.6) 

Hence using 

TrXX = —TrXX - TrX 2 , (2.7) 
dt 

we arrive at the following equation: 

^TrX 2 + /3^-TrX 2 = 0, (2.8) 
dt z dt 
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d o 
— TrX 2 
dt 



2TrX X , TrX 2 = TrX 2 . (2.9) 



t=o 

The solution of (2.8) is 

TrX 2 = 2TrX X ± ( 1 - e~ pt ) + TrX 2 . (2.10) 

It thus appears that whenever (3 > then the solution to (2.5) approaches the 
orbit such that 

TrX 2 = |TrX X + TrX 2 . (2.11) 
In the case with (3 = the solution (2.10) takes the form 

TrX 2 = 2TrX X t + TrX 2 , (2.12) 

i.e. there is no solution to (2.8) on the orbit unless 

TrX X = 0, (2.13) 

that is the time evolution starts from the orbit which is then by virtue of 
(2.12) and (2.13) the invariant set. Clearly, the same holds true for (3 < 0. 

Suppose that (3 > and TrX 2 > for all X G g. As we have shown above the 
solution to (2.5) approaches then the orbit whose parameters depend via (2.11) 
on the initial data. Furthermore, it is clear that (2.5) has the asymptotics 

TrX 2 

X + (3X + —X = (f x Ye' iX ) \^ x2=r2 - Y, (2.14) 
TrX 2 = r 2 , (2.15) 



where r = ^/|TrX Xo + TrX^. Evidently, (2.15) can be regarded as a con- 
straint. In section (4.2) (see also appendix) we demonstrate by an example of 
the SU(2) group that the solutions of the unconstrained equation (2.14) of 
the form 

TrX 2 

X + (3X + —X = (e«y e -«) | Trf2=r2 - y, (2.16) 

where r is a constant, approach the orbit TrX 2 = r 2 regardless of the initial 
data, that is this orbit is a universal attracting set. The only exception are 
the initial conditions such that 

X = iiY, X = vY, /x, v e R. (2.17) 

These initial data correspond to the following ansatz for the solution to (2.16): 

X(t) = x(t)Y, (2.18) 
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x + px H —xx 



Y 



'(e iX Ye~ iX ) 1 -y. (2.19) 

We note that (2.19) is not any asymptotics of the equation 

x 2 

x + f3x + — = 0, (2.20) 

x 

arising from inserting the ansatz (2.18) into (2.5). 

3 Periodic solutions to the Newton-like equation 

This section is devoted to the periodic solutions to the abstract Newton-like 
equation (2.5). Consider eq. (2.5). Our aim is to study the solutions to (2.5) 
on the orbit, therefore we write (2.5) as 

X + PX + ^-X = e iX Ye~ iX - Y, (3.1) 
c 

TrX 2 = const = c. (3.2) 
On using (3.1), (3.2) and the identity 

Tr(e iX Ye- iX - Yf = -2TrY '(e iX Y 'e~ iX - Y), (3.3) 
we arrive at the following relation: 

(TrX 2 ) 2 

TrX 2 + 2pTrXX + (3 2 TrX 2 - ± '- = -2TrY [e iX Y e~ iX - Y). (3.4) 

c 

Furthermore, an immediate consequence of (3.1) is 

TrX 2 

TrYX + pTrYX + = TrY{e iX Ye~ iX - Y). (3.5) 

c 

Assuming that 

= const = «, (3.6) 

we obtain from (3.5) 



c 

Inserting (3.7) into (3.4) we get 

2 ±) TrA" 2 - I 
c J c 



-TrX 2 = TrY[e lX Ye- lX - Y). (3.7) 



TrX 2 + 2/JTrXX + {(? + — ) TrX 2 - i(TrX 2 ) 2 = 0. (3.8) 
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TrY 2 7 \- -/ 

where F^O. Evidently, if (3.6) holds then 

V = X, (3.10) 

and (3.8) can be written as 

TrV 2 + 2pTrVV + f^ 2 + TrV 2 - ^-(TrV 2 ) 2 = 0. (3.11) 

Further, it follows immediately from (3.9) that 

TrVY = 0, TrV (n) Y = 0, (3.12) 
where is the n-th time derivative of V, and 

TrV 2 = TrX 2 - (T ^ )2 . (3.13) 
IrY z 

Clearly, whenever (3.2) and (3.6) take place then 

TrV 2 = const, TrVV = 0. (3.14) 

Suppose now (see (3.10)) that 

TrV 2 = TrX 2 = const = r, (3.15) 

and that the algebra g is three dimensional, then in view of (3.12), (3.14) and 
(3.15) the triple {Y, V, V} forms the orthogonal basis of g. On expanding V 
in this basis: 

V = fj,V + vY + pV, (3.16) 
and using a differential consequence of (3.15) such that 

TrVV = 0, (3.17) 

together with (3.12) we arrive at the following relation: 

V = pV (3.18) 

Hence using (3.14) one obtains 

TrVV = -TrV 2 = pTrV 2 . (3.19) 

Therefore, 

TrV 2 

V = -^V, (3^20) 
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TrV 2 

Finally, inserting (3.21) into (3.11) and taking into account (3.17) (see also 
(3.2), (3.6), (3.13) and (3.15)) we get 



1 



T 




o, 



(3.22) 



where e = TrY 2 . The validity of (3.22) is the condition for the existence of the 
solution to (2.5) satisfying (3.2), (3.6) and (3.15). In section 4.3 we show that 
in the case of the SU (2) group such the solution corresponds to the periodic 
motion on the sphere S2, more precisely, the uniform circular motion in a 
parallel. The counterpart of X is then the position vector in R 3 and (3.20) is 
simply the normal acceleration. Notice that the case r = corresponds to the 
equilibrium solutions to the equation (2.5) belonging to the orbit defined by 

(3.2) . In fact, taking into account (3.7) and (3.15) when r = as well as using 

(3.3) one finds 

X = XY. (3.23) 

Hence, by virtue of (3.2) we get 



A = ±J-. 



(3.24) 



The formula (3.23) related to the ansatz (2.19) is also obtained in the case 



c — 



which means that 



TrX 2 TrF 2 = (TiXY) 2 . 



(3.25) 



We end this section with the remark concerning equation (3.1), when /3 — 
and Y — 0. Clearly, we then have 

TrX 2 

where TrX 2 = const = c. Furthermore, an immediate consequence of (3.2) 
and (3.25) is 

TrX 2 = const. (3.27) 
In the light of the interpretation of the formula (3.21) in the case of the 
577(2) group mentioned above it is plausible that then (3.26) corresponds to 
the geodesic equations on the sphere S 2 - Obviously, the solutions to geodesic 
equations involve the periodic uniform motion in great circles. 



4 Dynamics on the SU(2) group 
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This section deals with the concrete realization of the abstract Newton-like 
equation (2.5) in the case of the SU(2) group. We first derive the nonlin- 
ear dynamical system generated by (2.5) satisfied by the group parameters. 
Consider eq. (2.5). On taking into account the following realization: 

J = §<r, (4.1) 

where cr = (<7i, a 2 , cr 3 ) and i — 1, 2, 3, are the Pauli matrices, of the 
infinitesimal generators Jj, i = 1, 2, 3, of the group SO(3) which is locally 
isomorphic to SU (2), we can write the general element of the Lie algebra X(t) 
as 

X(t)=x(t).|, (4.2) 
where x(t): R — > R 3 and the dot designates the inner product. Analogously 

Y = a-|, (4.3) 

where a is a constant vector of R 3 . Now, inserting (4.2) and (4.3) into (2.5) 
and using the elementary properties of the Pauli matrices we arrive at the 
following nonlinear system of second-order differential equations: 

x 2 sin |x| (a-x)x 



sin 


X 




X 





x + 3± H — -x = (cos Ixl — l)a H ;— r - ^ x x + (1 — cos Ixl)- — ^— , (4.4) 

x(0) = x , x(0) = x , 



where a x x designates the vector product of vectors a and x, |x| = vx 2 
stands for the norm of the vector x and x = ^. 

at 

Notice that the orbits given by (2.2) are the two-dimensional spheres 

x 2 = const. (4.5) 

The linear stability analysis shows that the equilibrium x = a, x = is 
unstable and the equilibrium x = —a, x = is stable for sufficiently large 
P > 0. 

We now discuss the asymptotic behaviour of the system (4.4). On making use 
of the relation 

TrAB = |a-b (4.6) 

for A = sl'y and B = b~, and taking into account (4.2) and (4.3) we arrive 
at the following realization of (2.10) in the case of the SU(2) group: 

x 2 = 2xo-xoi(l- e -^)+x 2 . (4.7) 
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2x -x + x^>0, (4.8) 
and P > 0, then the solution to (4.4) approaches the orbit 

x 2 = |x -x + x^. (4.9) 
Furthermore, it is clear that for the initial conditions such that 

xo-xo = 0, (4.10) 
and arbitrary (3 the orbit is an invariant set. 

The remaining solutions do not approach any orbit. Namely, for (3 > and 

|x -x + x 2 = 0, xq-xq^O, (4.11) 

the solutions to (4.4) tend asymptotically to the singular point x = 0. This 
point is approached after a finite period of time 

U = ~ln(l + -^V\, (4.12) 
p V 2x -x J 

for the initial data satisfying the inequality 

2x -x + x 2 <0. (4.13) 

Finally, if (5 < and x -x ^ then the trajectories go to infinity. 

Let us now focus our attention on the most interesting case when (3 > 
and the initial data satisfy the inequality (4.8) so that the solution to (4.4) 
approaches the sphere (4.9). We set for simplicity a = (0, 0, 03). Based on the 
observations of section 4.3 discussing asymptotic dynamics on the orbit (see 
formula (4.44)) we find that whenever the following condition holds: 

[3 2 V^xo-xo + x 2 , < \a 3 \ (l + cos y^xo-xo + x 2 ) , (4.14) 

where ft > and 

sin ^-bcp-XQ + x 2 
«3 — / ==- ^ 0, (4.15) 

^-ixo-xo + x 2 , 

then the system (4.4) has the limit cycle given by (see (4.43) and (4.45)) 



a 3 (l + cos \J 2(3 ^q-xq + x 2 , 
x j + x j = R ^ (4.17) 
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R 



(2/5- 1 x -x + x2) 



1 - 



/?4(2/5-lx -X +x2) 



+ xg) a 2 z 



\ (l + cos y^-ixo-xo 

On the other hand, if the following inequality is valid (see (4.48)): 

(3 2 ^2f3~ l x -± + xg > |a 3 | (l + cos y^-bto-xo + x§ ^ , 
then the trajectories tend to the equilibrium point such that 

xi = 0, x 2 = 0, ^ 3 = -sgn a 3 ^2p-hc^ic + x§, x = 0, 
where sgnx is the sign function. 



(4.18) 



(4.19) 



(4.20) 



It follows directly from (4.49) that the critical value of the bifurcation param- 
eter (3 is the implicit function of the initial data given by the equation 

/^/^xo-xo + xg = \a 3 \ (l + cos yfy^xo-xo + x§ ) ■ (4.21) 

It should be noted that in view of (4.14)-(4.18) the system (4.4) has an infinite 
number of limit cycles. Evidently, these limit cycles are unstable. 



4-2 Asymptotic dynamical system 

We now study the nonlinear dynamical system derived from (2.17) for the 
SU(2) group. In view of (4.4) such system is of the form 

sin T 1 — COS T 

x + P± + —x = (cosr - l)a H a x x H (a-x)x. (4.22) 

As with (4.4) we find that the equilibrium x = ±r x = 0, with the plus 
sign is unstable and the equilibrium with minus sign is stable for large enough 
(3>0. 

We now discuss the SU(2) realization of the equation (2.19). On setting 

x = ara, (4.23) 
we obtain from (4.22) the following equation: 



x + px H — -x x = (1 — cosr) —x — 1 . 



(4.24) 



10 



return to the asymptotic system (4.22). In appendix we show that whenever 
P > 0, then the solutions to (4.22) approach the orbit 



2 2 

x = r 



(4.25) 



for arbitrary initial data. The only exception are the initial conditions of the 
form 

x = /ia, x = ua, //, v e R, (4.26) 
corresponding to the case of the equation (4.24). 

Suppose that (3 > 0. As with the system (4.4) we set for simplicity a = 
(0,0, a 3 ). Suppose that the bifurcation parameter f3 satisfies the inequality 



(3< 



a 3 |(l + cosr) 



(4.27) 



where as^ 11 ^ 0. Using observations of section 4.3 we find that then the 
system (4.22) has the limit cycle of the form (see (4.43) and (4.45)) 



2^2 



^3 



(3 2 r 



a 3 (l + cosr) ' 
R\ 



1 - 



X ^ ~\~ X 2 

where the radius of the limit cycle (the circle (4.29)) is 

R = r, 

Further, for (3 such that 

(3> 



^ (1 + cos r) 2 a| 



(4.28) 
(4.29) 

(4.30) 



a 3 | (1 + cosr) 



(4.31) 



the solutions to (4.22) approach the equilibrium point 

X\ = 0, X2 = 0, x 3 = — sgna 3 r, x = 0. 
Evidently, the critical value of the bifurcation parameter (3 is given by 



(4.32) 



a 3 |(l + cosr) 



(4.33) 



From the discussion completed in section 4.3 it follows that the system (4.22) 
has the global Hopf bifurcation at /3 — f3 c . Phase portraits from numerical 
integration of the system (4.22) for (3 > are illustrated in Figs. 1 and 2. 
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4-3 Dynamics on the sphere S 2 



We now discuss the dynamics on the orbit, i.e. the two- dimensional sphere 
S2 which has been already mentioned above to be attracting or invariant set 
for the investigated nonlinear dynamical systems for (5 > 0. In view of (4.22) 
the corresponding system can be written as 

IT^ Sill T 1 — COS T 

ii + (3ii H — -n = (cosr — l)a H a x n H (a-n)n, (4.34a) 

n 2 = r 2 . (4.34b) 

Clearly eq. (4.34b) is a constraint on (4.34a). By switching over to the spher- 
ical coordinates such that 

n = (r sin 9 cos tp, r sin 9 sin ip, r cos 9) , (4.35) 

and setting (without loose of generality) 

a= (0,0, a 3 ), (4.36) 

we obtain from (4.34) the following system: 

... ]^ — COS T 

9 + /39-sm9 cos 9 p 2 - a 3 sin 9 = 0, (4.37a) 

r 

sin t 

(p + p<p+ 2ctg^ 9(p - a 3 = 0. (4.37b) 

r 

We note that the system (4.37) is Lagrangian one. The Lagrangian is 

L = \{9 2 + sin 2 9 p 2 ) - a 3 l ~ COST cos 9, (4.38) 
2 r 

and the system (4.37) can be written in the form of the Lagrange equations 

ddL dL 

d dL dL ^ . t . 

Its* - T v = «*• < 439) 

where the generalized nonpotential forces are 

sin t 

Q e = -(39, Q v = -13 sin 2 9 <p + a 3 sin 2 ^. (4.40) 
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(4.39) as a generalized spherical pendulum. The standard spherical pendulum 
is easily seen to correspond to the particular case of /3 — and r = tt. The 
case (5 — and r = 2ti or j3 = and a 3 = refers to the free motion on the 
sphere. 

Consider the system (4.37). Let us assume that a^^ 1 ^ 1 - ^ 0. On setting 9 = 
const, we find tp = uo = const, so we then deal with the uniform circular 
motion in a parallel. Furthermore, we have 

2 1 COS T . A _ . 

- u 2 = a 3 — , 4.41a 

r cos 9 

fa = a 3 — . (4.41b) 

r 

By virtue of (4.41a) 9 G (f,7r] ("South of the equator") for a 3 > 0, and 
9 G [0, |) ("North of the equator") for 03 < 0. The singular case 9 — | can be 
treated as the motion in the equator with infinite angular velocity. In view of 
(4.41b) such a case corresponds to the limit (3 = 0. 

We remark that the abstract formula (3.22) providing the criterion of periodic 
motion on the orbit takes the form 

uo 2 cos 2 9 + f? + = 0, (4.42a) 

r 

where uj — ip — const and 9 = const. On the other hand, the abstract relations 
(3.7) and (3.6) taken together yield 

, 1 — cosr , 4 . 

u = -°3 7T> 4 - 42b 

r cost^ 

i.e. the formula (4.41a) is obtained. It can be easily checked that the system 
(4.41) is equivalent to the system (4.42). 

We now return to the system (4.41). Let us assume that (3 7^ and r tt, 
so we exclude the case of the standard spherical pendulum. An immediate 
consequence of (4.41) is the formula 

3 2 r 

cos 9 = -. (4.43) 

a 3 (l + cosr) 

Hence, recalling that the case of (3 > is investigated, we find 



P < J |Q3|(1 + COSr) . (4.44) 
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We now come to the analysis of the asymptotic dynamics on the sphere S 2 . 
Consider the system (4.37). Integrating (4.37b) yields 




— / sin 2 9 e^dr + y? sin 2 9 ) 4^ . (4.46) 
r J / sin 9 

o / 



We discuss the most interesting case of f3 > 0. Assuming that sin6> y4 we 
obtain from (4.46) the following asymptotic relation: 

/ sin 2 9 e^ T dr 
smr 5 smr 1 

V? = a 3 hm - — ; — 7t— — -T- = a 3 ■ . (4.47) 

r *->°o sin 2 e# r 2ctg# 9 + (3 

Suppose now that g^^ 1 ^ 0. Inserting (4.47) into (4.37b) we find that at 
asymptotic times <p = 0, i.e. (p = const, which in view of (4.47) leads to 9 = 
const. Clearly, the case sin 9 — > corresponds to approaching the equilibrium 
point 9 = it for a 3 > or 9 = for a 3 < 0. Indeed, the equilibrium 9 = 
for a 3 > and 9 = it for a 3 < have been already mentioned in section 
4.1 to be unstable. Taking into account (4.41) we see that whenever (3 > 
and (13^- 7^ then for (3 satisfying (4.44) the trajectories approach the circle 
given by (4.43). On the other hand, whenever (3 fulfils the inequality 



(3 > J Q3|(1 + C ° Sr) , (4.48) 



then the trajectories tend to the equilibrium point 9 = n for a 3 > and 9 = 
for a 3 < 0. Therefore, the critical value of the bifurcation parameter (3 is 



ft = ^M!±£5!I>. (449) 

Notice that for (3 = (3 C the radius of the circle given by (4.45) equals zero, 
that is the periodic trajectory (limit cycle) reduces to the equilibrium point 
9 = ix {9 = 0). Furthermore, as we have shown above the limit cycle given 
by (4.43) and the equilibrium point 9 = n {9 = 0) are universal attracting 
sets for < (3 < (3 C and (3 > f3 C) respectively. We conclude that whenever 
a3 smr q an j p > q ^j ien th e system (4.37) has a global Hopf bifurcation 
(the bifurcation of a limit cycle from an equilibrium point) at (3 — (3 C . 
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find that the system (4.37) is asymptotically Hamiltonian one. 



5 Conclusion 



In this work we have introduced an abstract second-order Newton-like dif- 
ferential equation on a general Lie algebra such that orbits of the Lie-group 
action are attracting set. By passing to the group parameters corresponding 
to the concrete Lie algebra the Newton-like equation generates the nonlinear 
dynamical system referred to as the group-invariant one with the attracting 
set coinciding with the orbit. In other words, a method has been found in 
this paper for constructing nonlinear systems with the presumed form of an 
attracting set. We have demonstrated that in the case of the SU(2) group the 
attracting manifold is the universal attracting set for the system generated by 
the asymptotic version of the abstract Newton-like equation. It is suggested 
that such universality holds true, at least in some regions of a phase space, for 
general Lie groups. In order to clarify the asymptotic behaviour of the non- 
linear dynamical systems arising in the case with the SU{2) group we have 
studied the dynamical system on the sphere S 2 corresponding to the general- 
ized spherical pendulum. Such system has been shown to have the global Hopf 
bifurcation. Based on this observation we have found that the group-invariant 
dynamical system is an interesting example of the system which has an infinite 
number of unstable limit cycles. On the other hand, it has turned out that the 
asymptotic dynamical system has a global Hopf bifurcation. We note that the 
discussed limit cycles (circles) arising in the case of the SU (2) group coincide 
with sections of the sphere S 2 by the planes perpendicular to the vector a. In 
our opinion, such a coincidence takes place in the general case of an orientable 
two-dimensional Lie-group manifold, that is the concrete form of the limit cy- 
cles (whenever they exist) is determined by the orbit of the Lie-group action 
and the direction of the element Y . Thus we suggest that the actual treat- 
ment provide a method for constructing nonlinear dynamical systems with the 
presumed form of a limit cycle. The example of the SU (2) group discussed 
herein with simple oscillatory dynamics on the orbit is a model one. We re- 
call that numerous group-theoretic constructions are tested by the example of 
the SU(2) group. It seems that in the case of groups with higher-dimensional 
orbits the dynamics shoul be much more complex. Since we do not know any 
criterion excluding symmetry of the chaotic attractor, therefore we should ad- 
mit that the dynamics would be chaotic one. We also observe that even in the 
case with the two-dimensional orbits some open questions arise when the orbit 
is not a connected set. For example, it is not clear referring to the experience 
with the SU(2) group what is the behaviour of the system having the orbit 
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of the above comments it seems that the group-theoretic approach introduced 
herein would be a useful tool in the theory of nonlinear dynamical systems, 
especially in the bifurcation theory, as well as in such branches of applied 
mathematics as for example the control theory. 
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We now demonstrate that the solutions to the system (4.22) with f3 > ap- 
proach the sphere x 2 = r 2 for arbitrary initial data except those ones satisfying 
the condition (4.26) referring to the case of the equation (4.24). Consider the 
system (4.22). We assume for simplicity (without loose of generality) that 
a = (0, 0, as). By switching over to the spherical coordinates such that 

Xi = p sin # cosy?, 
X2 = p sin 9 sin <p, 

x 3 = pcos9, (A.l) 
we arrive at the following system: 

p + (3p + p(9 2 + sin 2 9 if 2 ) - ?j + p 2 ^ = a 3 (l - cosr) cos 9 ^ -(A^a) 

2,0 . 1 cos v 

9 + {39-sm9 cos 9 y? 2 + — 9 - a 3 sin 9 = 0, (A.2b) 

P P 
2 sm v 

(3 + f3ip + 2ctg9 9(p + —(p-a 3 = 0. (A.2c) 

p r 

Integration of (A. 2c) yields 

[as 8 -^ I P 2 sin 2 9 e^dr + ip p 2 sin 2 O ) 

V = ± — — r^r- a '- • (A.3) 

p 2 sm 2 9 

Suppose that (3 > and sin9 y4 0. We also assume that as^ 1 - ^ 0. Notice 
that this inequality implies p -/-> since then x = is not a solution of 
(4.22). Proceeding analogously as with eq. (4.46) we arrive at the following 
asymptotic relation: 

Sil1 r 1 I A A \ 

m = a 3 ^ = . (A.4) 

r J + 2ctg99 + (3 

In view of (A. 2c) this leads to (p — 0, i.e. ip = const. Hence we find that at 
asymptotic times 

- + ctg9 9 = const. (A.5) 
P 

Integrating (A.5) we obtain 

psm9 = Ce^, (A.6) 

where C and p are constant. Suppose that p > 0. By virtue of (A.6) this 
implies p — > oo. Hence, taking into account (A.2a)-(A.2c) and differential 
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p z tg z 9{pp-f>y p z sm6<p 2 1-cosr 2 
- /fyi r r h a 3 5 P cos 9 = p . (A. 7) 

But p2sit f/^ 2 > a 3 i^^p cos 9. Therefore, the left-hand side of (A.7) is neg- 
ative and the right-hand side is positive. This contradiction shows that p in 
(A. 6) cannot be positive. Let us assume now that p < 0. Then p — > and in 
view of (A.l) x -> 0. But x = is not the solution to (4.22) when a 3 ^ ^ 0. 
Thus it turns out that p — and the asymptotic formula (A. 6) reduces to 

psin9 = C. (A.8) 

Now, eqs. (A.2a)-(A.2c) and (A.8) taken together yield 

- o P ^l n + (r 2 - p 2 sin 2 eW + a 3 (l - cos r)p cos # = 0. (A.9) 
p 2 cos 2 

On integrating (A.9) written with the help of (A.8) in the form independent of 
9, one can easily find that p ^ at asymptotic times, leads to contradiction. 
Therefore p = 0, i.e. p = const. Hence by virtue of (A.8) we have 9 = const. 
Assuming that p ^ r, using (A. 2a) and (A. 2b), and taking into account that 
p 7^ 0, sin 9 7^ 0, (p 2 ^ 0, we get pp 2 = 0, a contradiction. Thus we find p = r. 
Notice that an immediate consequence of this observation and eq. (A.9) is the 
formula (4.41a) describing the uniform circular motion on the sphere x 2 = r 2 . 
Clearly, the case sin 9 — > corresponds to approaching the stable equilibrium 
x = (0, 0, — sgna 3 r), x = of (4.22). We finally note that in the case of the 
ansatz (4.23) implying eq. (4.24) we have sin# = in (A.l) and the system 
(A. 2) becomes singular one. 
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Fig. 1. The system (4.22) with (5 = 0.1, r = 1, a = (0, 0, 1) and x = (1.4, 1, 1), 
xo = (0.1,0.1,0.1). The parameter f3 satisfies (4.27). The initial condition 
fulfils Xq > r 2 . The trajectory starts from the point marked with the cross. 
Left: the projection on the (xi,x 2 ) plane. Right: the projection on the (xi,x 3 ) 
plane. 

Fig. 2. The system (4.22), where /3, r, a, are the same as in Fig. 1, and 
x = (0.1,0.1,0.1), x = (0.5,0.5,0.5). The initial data satisfy x 2 , < r 2 . Left: 
the projection on the (x±, x 2 ) plane. Right: the projection on the (x±, x 3 ) plane. 



20 



